HAAC 

Joumals.ASM.org 

Population Pharmacokinetic Assessment of the Effect of Food on 
Piperaquine BioavailabiUty in Patients with UncompUcated Malaria 

Joel Taming,'''' Niklas Lindegardh,''''t Khin Maung Lwin," Anna Annerberg,' '' Lily Kiricharoen," Elizabeth Ashley,^ ''''' 
Nicholas J. White,' '' Francois Nosten,' '''^ Nicholas P. J. Day' '' 

Mahidol Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University, Bangl<ok, Thailand'; Centre for Tropical Medicine, Nuffield Department 
of Medicine, University of Oxford, Oxford, United Kingdom''; Shoklo Malaria Research Unit, Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, 
Mahidol University, Mae Sot, Thailand^ 

Previously published literature reports various impacts of food on the oral bioavailability of piperaquine. The aim of this study 
was to use a population modeling approach to investigate the impact of concomitant intake of a small amount of food on piper- 
aquine pharmacokinetics. This was an open, randomized comparison of piperaquine pharmacokinetics when administered as a 
fixed oral formulation once daily for 3 days with (n — 15) and without {n — 15) concomitant food to patients with uncompU- 
cated Plasmodium falciparum malaria in Thailand. Nonlinear mixed-effects modeling was used to characterize the pharmacoki- 
netics of piperaquine and the influence of concomitant food intake. A modified Monte Carlo mapped power approach was ap- 
plied to evaluate the relationship between statistical power and various degrees of covariate effect sizes of the given study design. 
Piperaquine population pharmacokinetics were described well in fasting and fed patients by a three-compartment distribution 
model with flexible absorption. The final model showed a 25% increase in relative bioavailability per dose occasion during re- 
covery from malaria but demonstrated no clinical impact of concomitant intake of a low-fat meal. Body weight and age were 
both significant covariates in the final model. The novel power approach concluded that the study was adequately powered to 
detect a food effect of at least 35%. This modified Monte Carlo mapped power approach may be a useful tool for evaluating the 
power to detect true covariate effects in mixed-effects modeling and a given study design. A small amount of food does not affect 
piperaquine absorption significantly in acute malaria. 



Malaria is the most important parasitic disease in humans and 
kiUs approximately 660,000 people each year (1). Artemis- 
inin-based combination therapy (ACT) is the recommended first- 
line treatment worldwide. The artemisinin component is a very 
potent but short-acting drug and kills the majority of parasites, 
while the longer-acting partner drug is slowly eliminated and kiUs 
residual parasites to prevent recrudescent malaria. 

The fixed oral combination of dihydroartemisinin and piper- 
aquine is a promising ACT, which has shown excellent cure rates 
of 98.7% (95% confidence interval [CI], 97.6 to 99.8%) in a 
pooled analysis of 3,547 patients with malaria from 12 different 
sites (2). Dihydroartemisinin-piperaquine is currently recom- 
mended as a once-daily dose for 3 days based on patient body 
weight (daily target dose of 18 mg of piperaquine phosphate/kg of 
body weight) (3). However, concerns have been raised that small 
children are underdosed due to the nonlinear relationship be- 
tween elimination clearance and body weight (4). 

Dihydroartemisinin has a very short terminal elimination half- 
life of approximately 1 h, while piperaquine has a long half-life of 
approximately 20 to 30 days (4-7). The absolute oral bioavailabil- 
ity of piperaquine has not been reported in humans, since paren- 
teral formulations are unavailable. However, a 50% oral bioavail- 
ability compared to an experimental parenteral formulation was 
reported in rats ( 8 ) . Piperaquine is a highly lipophilic drug, and its 
absorption might be facilitated by concomitant intake of fat, as 
described previously for other lipid-soluble antimalarial drugs 
such as lumefantrine (9), halofantrine (10), mefloquine (11), and 
atovaquone (12). 

Contradictory results have been reported for the effects of con- 
comitant food intake on piperaquine pharmacokinetics in healthy 
volunteers and patients with malaria. Sim and colleagues reported 



a 98% increase in total exposure of piperaquine after a high-fat 
breakfast (53 g fat) in healthy Caucasian volunteers (;i = 8; cross- 
over design) (13). Nguyen et al. reported a more modest 41% 
increase in total piperaquine exposure after a standardized Viet- 
namese breakfast (17 g fat) in healthy Vietnamese volunteers (« = 
14; parallel design) (14). However, Hai and colleagues reported no 
significant impact on piperaquine pharmacokinetics after a simi- 
lar standardized Vietnamese breakfast (17 g fat) in healthy Viet- 
namese volunteers (« = 32; parallel design) (15). The noncom- 
partmental analysis of data presented here showed no effect of 
concomitant intake of a small amount of fat (6.4 g fat) on piper- 
aquine exposure in patients with uncomplicated Plasmodium fal- 
ciparum malaria in Thailand (n = 30; parallel design) (16). How- 
ever, all reported studies used a noncompartmental approach 
with a low statistical power for detecting concomitant food intake 
as an influential factor on piperaquine pharmacokinetics, com- 
pared to a modeling approach. 

The aim of this study was to use a potentially more powerful 
population approach to analyze data from a previously reported 
study (16) to evaluate the effects of concomitant food intake on 
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TABLE 1 Demographics of patients with uncomphcated P. falciparum 
malaria in Thailand 



Parameter 


Value for group 




Fasting patients 


Fed patients 


Pharmacokinetic data 






Total no. of patients 


15 


15 


Total no. of piperaquine samples 


535 


541 


Median daily dose of piperaquine 


17.2 (16.0-18.6) 


17.5 (16.0-18.8) 


phosphate (nig/kg) (range) 






Median daily dose of 


2.14 (2.00-2.33) 


2.19 (2.00-2.35) 


dihydroartemisinin (mg/kg) 






(range) 






Demographics 






Median age (yi') (range) 


38 (18-55) 


28 (19-45) 


Median body wt (kg) (range) 


50 (39-62) 


53 (45-73) 


No. of males/no. of females 


13/2 


13/2 


Median axillary temp at 


36.5 (36.2-38.2) 


37.1 (35.9-39.2) 


admission (°C) (range) 






Median parasitemia at admission 


8,000 (448-140,000) 


8,000 (352-60,000) 


(no. of parasites/|xl) (range) 






Median diastolic blood pressure 


70 (60-80) 


70 (60-110) 


(mmtig) (range) 






Median systolic blood pressure 


110 (90-130) 


110 (90-140) 


(mmHg) (range) 






Median hematocrit (%) (range) 


41 (30-45) 


42 (33-47) 


Median pulse (beats/min) (range) 


80 (65-96) 


84 (72-120) 



the pharmacokinetic properties of piperaquine in patients with 
uncomphcated P. falciparum malaria in Thailand. 

MATERIALS AND METHODS 

Study design. This was an open, randomized, parallel study of piper- 
aquine pharmacokinetics when administered as a fixed oral formulation 
once daily for 3 days with and without concomitant fat to patients with 
uncomplicated P. falciparum malaria in Thailand. Clinical and noncom- 
partmental results were reported in full previously (16). Briefly, 30 pa- 
tients aged 16 to 65 years were enrolled and randomized into one of the 
two treatment arms. Inclusion criteria were microscopy confirmation of 
asexual P. falciparum or mixed infections, no signs of severe malaria, and 
willingness to participate in the study (fully informed consent). Exclusion 
criteria were a4% red blood cell parasitemia, positive urine test for preg- 
nancy, known hypersensitivity to artemisinins or piperaquine, treatment 
with dihydroartemisinin or piperaquine within the past 4 months, or 
hematocrit level of <30%. Full demographics are given in Table 1. Ethical 
approval was granted by the Faculty of Tropical Medicine Mahidol Uni- 
versity Ethical Committee, Bangkok, Thailand, and the Oxford Tropical 
Research Ethics Committee (OxTREC), United Kingdom. 

Drug admmistration. All patients were given a supervised standard 
fixed oral formulation of dihydroartemisinin and piperaquine (40 mg 
dihydroartemisinin and 320 mg piperaquine phosphate per tablet) (Duo- 
Cotecxin; Beijing HoUey-Cotec Pharmaceuticals Co., Ltd., China). A 
weight-based dose regimen, once-daily treatment for 3 days, was em- 
ployed to achieve a daily target dose oflSmgof piperaquine phosphate/kg 
ofbody weight (Table 1). All patients in the fed group (n = 15) received a 
200-ml carton of chocolate milk, containing 6.4 g of fat, with each dose. 
Patients in the fasting group (n = 15) were given study medication at 
enrollment, after at least 2 h of fasting, and consecutive doses after an 
overnight fast. Patients were asked to continue fasting for 3 h after each 
dose. 

Pharmacokinetic sampling. Blood samples during the intensive col- 
lection phase (up to 12 h after the last dose) were collected through an 
indwelling intravenous cannula flushed with 0.5 ml heparinized saline 
solution after each sample collection. A total of 0.5 ml of blood was dis- 
carded immediately before sampling, and blood samples (2 to 5 ml) were 
collected into lithium heparin tubes predose (0 h); at 0.5, 1, 2, 3, 4, 7, and 



24 h after the first dose; at 1, 3, 4, 5, 7, and 24 h after the second dose; and 
at 1, 2, 3, 4, 5, 6, 8, and 12 h after the third dose. Additional samples were 
collected by venipuncture on days 4, 5, 7, 14, 2 1 , 28, 42, 56, 70, 84, 98, 112, 
and 126 after the first dose. All blood samples were centrifuged at 2,000 X 
g for 10 min, and plasma was stored in liquid nitrogen within 30 min of 
collection. Plasma samples were transported on dry ice to the Department 
of Clinical Pharmacology at the Mahidol-Oxford Tropical Medicine Unit, 
Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand. 
Piperaquine plasma samples were quantified by using a previously re- 
ported high-throughput assay consisting of liquid chromatography linked 
with tandem mass spectrometry detection (17). Triplicates of quality con- 
trol samples at three different concentrations (i.e., 4.5, 20, and 400 ng/ml) 
were analyzed within each batch of samples to ensure precision and accu- 
racy during drug measurements. The coefficient of variation was below 
5% for all quality control samples. The lower limit of quantification was 
1 .2 ng/ml. 

Population pharmacokinetics. Piperaquine plasma concentrations 
were transformed into their natural logarithms and evaluated by using 
nonlinear mixed-effects modeling in NONMEM v.7.2 (Icon Develop- 
ment Solutions, MD). Automation and visualization of diagnostic results 
were performed by using Pearl-Speaks-NONMEM (PsN) v.3.6.2 (18, 19), 
Xpose V.4 (20), R v.2.13.1 (R Foundation for Statistical Computing), and 
Piranav.2.6 (21). 

Subroutine ADVAN5 and the first-order conditional estimation 
method with interactions were used throughout modeling. The objective 
function value (OFV) computed by NONMEM as minus twice the log 
likelihood of data, goodness-of-fit diagnostics, and simulation-based di- 
agnostics were used to discriminate between models. A total of 91 out of 
1,016 samples (9.1%) were quantified below the limit of quantification 
and omitted during model building. The majority of these omitted sam- 
ples (82.6%) were in the late terminal phase (>70 days) during the long 
foUow-up period. The impact of omitting these samples was evaluated by 
using the final model. 

One-, two-, three-, and four-compartment distribution models were 
first evaluated during model building. The best-performing distribution 
model was assessed together with different absorption models (i.e., first- 
order absorption with and without lag time and a more flexible transit 
compartment absorption with a fixed number of transit compartments 
for the population [22] ). Implementation of a fixed relative bioavailability 
of 100% for the population but allowing for between-subject variability 
and between-dose occasion variability (i.e., modeled as within-subject 
variability between doses) in the same parameter was also evaluated. Phar- 
macokinetic parameters were modeled assuming log-normal distribu- 
tions, while between-subject variability and between-dose occasion vari- 
ability were assumed to be normally distributed with a zero mean and 
variance. The random residual variability was assumed to be additive, 
since data were modeled as natural logarithms (i.e., essentially equivalent 
to an exponential error model for untransformed data). 

Body weight was investigated as a covariate by using an aUometric 
function with power values of 3/4 for clearance parameters and 1 for 
volume parameters. Systematic variability in the relative bioavailability 
between dose occasions was evaluated as a categorical and a linear cova- 
riate for doses 1, 2, and 3. Effects of age, initial parasitemia, hematocrit, 
temperature upon admission, sex, and concomitant food intake on all 
parameters were evaluated formally with a stepwise forward-addition 
(P< 0.05) and backward-deletion (P< 0.01) covariate approach by using 
the automated Stepwise Covariate Model (SCM) implemented in PsN. 
Continuous covariates were tried as linear, piecewise linear, power, and 
exponential functions, and categorical covariates were tried as linear func- 
tions. The final structural model, including body weight as an aUometric 
function on clearance and volume parameters as well as between-subject 
variability on all parameters and between-dose occasion variability on 
relative bioavailability and mean absorption time, was also used to inves- 
tigate the impact of concomitant fat intake on absorption parameters with 
a fuU covariate approach (23). Concomitant fat intake was simultaneously 
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implemented as a categorical food effect on the relative bioavailability and 
absorption parameters, followed by bootstrap diagnostics {n = 500) to 
evaluate the clinical relevance of such potential covariate effects. 

Basic goodness-of-fit diagnostics, visual and numerical predictive 
checks, as well as bootstrap diagnostics were used to evaluate the appro- 
priateness of the final model. Visual predictive checks were visualized by 
plotting the 95% confidence intervals of the 5th, 50th, and 95th simulated 
(n = 2,000) percentiles against measured piperaquine concentrations. 
The nonparametric bootstrap diagnostics {n = 1,000) were stratified by 
drug administration with or without food to maintain an equal distribu- 
tion of patients in the resampled data. 

Evaluation of population versus noncompartmental approaches. A 
traditional power calculation was performed by using the mean results ( ± 
standard deviations [SD]) of total piperaquine exposure from the non- 
compartmental analysis reported previously (16). The mean total piper- 
aquine exposure in fasting patients was used to calculate the mean expo- 
sure in fed patients after a putative food effect, and the SD was based on 
that observed for fasting patients but was assumed to be proportional to 
exposure. Standard functionalities in STATA v. 12.0 (Stata Corp., TX, 
USA) were used to calculate the minimum effect size of food needed for a 
statistically significant difference between a fasting group {n = 15) and a 
fed group {n = 15) (power = 0.8; alpha = 0.05). 

An embedded functionality in PsN for rapid sample size calculations 
of mixed-effects models (24) was modified to calculate the minimum 
effect size of food on relative bioavailability needed for a statistically sig- 
nificant covariate relationship (modified Monte Carlo mapped power ap- 
proach). The original study design of 15 fed and 15 fasting patients was 
extended to include 1 1,000 simulated patients (50% in the fed group and 
50% in the fasting group). Individual concentration-time data were sim- 
ulated for these patients with the final model with the addition of a fixed 
categorical effect size of 5% to 50% of food on relative bioavailability (i.e., 
1,000 patients for each effect size with an increment of 5% between sim- 
ulations). The simulated data were then separately reestimated for each 
effect size with the final model with and without a food effect on relative 
bioavailability. The difference in individual objective function values 
(AiOFV) was calculated for each patient for the two models. Fifteen 
AiOFVs were resampled at random (bootstrap; n = 10,000) from each 
group (i.e., fed and fasting) and effect size (i.e., 5% to 50%), and fractions 
of the 10,000 bootstrap samples above the critical total OFV (2 AiOFV) of 
3.84 and 6.63 (i.e., power at alpha values of 0.05 and 0.01, respectively) 
were plotted against the evaluated covariate effect sizes. 

RESULTS 

The study medication was well tolerated, with no severe adverse 
reactions reported during the study period. There were no statis- 
tical differences between the two groups in terms of demographic 
variables (Table 1). A total of eight patients had recurrent malaria 
during the 126 days of follow-up, but none of these patients were 
classified as having recrudescent malaria by PGR genotyping. The 
median time to a new infection was 42 days (range, 32 to 87 days) 
in the fasting group (n = 5) and 70 days (range, 58 to 125 days) in 
the fed group (n = 3). 

Population pharmacokinetics. Piperaquine population phar- 
macokinetics were well described by a 3-compartment distribu- 
tion model, with no significant improvement by adding an addi- 
tional peripheral distribution compartment (AOFV = —5.08). A 
transit compartment (n = 3) absorption model described the ab- 
sorption phase well and was superior to all other absorption mod- 
els (AOFV < —51.9). The absorption rate from the last transit 
compartment could be set as identical to the rate constant between 
transit compartments without a significant impact on the model 
(AOFV = 0.41). This also provided a more stable absorption 
model. The final structural model is shown in Fig. 1. 

Implementation of a fixed relative bioavailability of 100% for 



n = 3 



Gut 
compartment 
(F X dose) 



r 




ktr 










( s 

Transit (3) 
compartment 



Periplieral 
Compartment 



Q1/V0 



ktr 



Centrai 
Compartment 
(Vc/F) 



CL/Vr 



Q2/ Vp. 



Q2/VC 



Periplieral 
Compartment 



FIG 1 Final structural model for piperaquine population pharmacokinetics in 
fasting (« = 15) and fed (n = 15) patients with uncomplicated P. falciparum 
malaria in Thailand. CL, elimination clearance; F, relative oral bioavailability; 
ktr, transit absorption rate constant; Q, intercompartment clearances; V^, ap- 
parent volume of distribution of the central compartment; Vp, apparent vol- 
ume of distribution of the peripheral compartments; MTT, mean absorption 
transit time; n, number of transit compartments [MTT = (n + l)lktr]. 



the population but allowing for between-subject variability and 
between-dose occasion variability in the same parameter im- 
proved the model fit significantly (AOFVs of —46.1 and —160, 
respectively, when implemented sequentially). However, the be- 
tween-subject variability could be removed after implementation 
of between-dose occasion variability without a major impact on 
the model fit (AOFV = 1.94). Between-dose occasion variability 
in the mean absorption time also improved the model signifi- 
cantly (AOFV = -27.6). 

A fixed allometric function for body weight improved the 
model fit marginally (AOFV = —3.68) but was kept in the final 
model based on prior strong physiological evidence for such a 
covariate relationship. Dose-occasion as a categorical covariate 
effect on relative bioavailability improved the model significantly 
(AOFV = -11.0) and resulted in 21.9% and 50.9% increased 
bioavailability at dose 2 and dose 3, respectively, compared to dose 
1. However, this could be simplified to a linear covariate relation- 
ship (i.e., 25.3% increase in bioavailability per dose) without a 
significant reduction in model fit (AOFV = —0.062). Effects of 
age and sex on peripheral distribution volume and concomitant 
food intake on mean absorption time were significant covariates 
in the forward-addition step {P < 0.05), but only the effect of age 
on the peripheral volume of distribution could be retained in the 
backward-elimination step with a more parsimonious cutoff {P < 
0.01). This covariate relationship resulted in a linear increase in 
the peripheral volume of distribution of 4.10% for each year of age 
increase. A separate bootstrap diagnostic from the implementa- 
tion of a food effect on mean absorption time and relative bio- 
availability showed a significantly higher absorption rate during con- 
comitant food intake but no effect on total bioavailability (Fig. 2). 

Basic goodness-of-fit diagnostics resulted in adequate model 
performance, with no obvious model misspecification (Fig. 3). 
However, a minor trend of data censoring at the lower limit of 
quantification could be seen, but this resulted in no model mis- 
specification in terms of simulated and observed fractions of cen- 
sored data (data not shown). Furthermore, omitting data below 
the limit of quantification resulted in almost identical parameter 
estimates as when these were handled as categorical data (i.e., the 
M3 method resulted in a <5.5% absolute mean bias compared to 
when the data were omitted) (25). The predictive checks demon- 
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FIG 2 Box plots (interquartile ranges with 2.5 to 97.5 percentiles) showing the 
effect of estimated food effects on mean absorption transit time (MTT) and rela- 
tive bioavailability (F). Vertical dashed lines indicate no effect and ±25% effects. 



strated excellent predictive performance of the final model (Fig. 4) 
and calculated 4.32% (95% CI, 2.38 to 8.43%) and 5.62% (95% 
CI, 2.27 to 8.32%) of the observed data below and above the sim- 
ulated 90% prediction interval. Final parameter estimates and 
bootstrap diagnostics are summarized in Tables 2 and 3. 



Evaluation of population versus noncompartmental ap- 
proaches. A noncompartmental analysis and a groupwise statisti- 
cal comparison of total piperaquine exposure between the fed and 
fasting groups resulted in 80% power (alpha = 0.05) to detect a 
minimum difference of 81%. However, the novel power method- 
ology for mixed-effects models resulted in 80% power (alpha = 
0.05) for a covariate effect of 35%, thus proving the study design to 
be adequately powered for a population approach if a putative 
clinically relevant difference of >35% would be present between 
fed and fasting patients. The statistical power to detect various 
degrees of covariate effects is illustrated in Fig. 5. 

DISCUSSION 

The fixed oral combination of dihydroartemisinin and piper- 
aquine has demonstrated excellent cure rates in adult patients 
with P. falciparum malaria in Thailand (2). In this small study, 
there was no recrudescent malaria when the combination was ad- 
ministered alone or together with a standardized meal. The rela- 
tively high rate of new infections is most likely a consequence of 
the long follow-up period of 125 days. 

Population pharmacokinetics. The population pharmacoki- 
netics of piperaquine were well characterized by the study design 
and the developed nonlinear mixed-effects model. Piperaquine 
showed a multiphasic distribution, as described previously (4, 7, 
26). A more flexible absorption model improved the model sig- 
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FIG 3 Goodness-of-fit diagnostics of the final population pharmacokinetic model of piperaquine in fasting (fi = 15) and fed (« = 15) patients with uncom- 
plicated P. falciparum malaria. Broken lines, locally weighted least-squares regression; solid lines, line of identity; broken horizontal lines, lower limit of 
quantification. The observed concentrations, population predictions, and individual predictions were transformed into their logarithms (base 10). 
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FIG 4 Visual predictive clieck of the final model describing the population 
pharmacokinetics of piperaquine in fasting (n = 15) and fed (« = 15) 
patients with uncomplicated P. falciparum malaria. The inset shows piper- 
aquine simulations at 0 to 72 h. Open circles, observed data points; solid lines, 
5th, 50th, and 95th percentiles of tlie observed data; shaded areas, 95% confi- 
dence intervals of simulated (« = 2,000) 5th, 50th, and 95th percentiles. Bro- 
ken horizontal lines are the lower limit of quantification. Venous plasma pip- 
eraquine concentrations were transformed into their logarithms (base 10). 

nificantly and described the variable absorption of piperaquine 
satisfactorily. A major improvement was seen when allowing for 
between-dose occasion variability, which confirms the large vari- 
ability in the absorption of piperaquine both between and within 



subjects. The estimated relative bioavailability of piperaquine in- 
creased by 25.3% for each consecutive dose. This was not related 
to concomitant food intake and might be a consequence of pa- 
tients recovering from malaria during the course of the treatment 
and an improvement in general clinical status. This phenomenon 
has been described previously for piperaquine treatment in preg- 
nant women (7) as well as for mefloquine treatment in children 
with acute uncomplicated P. /aZapanim malaria (27). The periph- 
eral volume of distribution increased with age in this study, which 
is likely to reflect the increased body fat-to-water content in older 
patients. However, this is unlikely to have clinical implications 
considering that total exposure of piperaquine was not affected. 
Final pharmacokinetic parameter estimates were generally in 
good agreement with values reported for adult patients with ma- 
laria (5, 7). 

This study demonstrated no clinical impact on concomitant 
intake of a small amount of fat on the total absorption, day 7 levels, 
peak concentrations, or time to peak concentrations of piper- 
aquine. Previous studies in healthy volunteers reported a large 
impact on the relative bioavailability of piperaquine after a high- 
fat meal (13,14). These differences are likely to be explained by the 
small amount of fat administered in this study, but we believe that 
it reflects the clinical reality, as many patients are unable to con- 
sume a fall meal when presenting with malaria. Most patients are 
unwell when presenting with malaria, and nausea is commonly 
seen in these patients. One carton of milk seems a reasonable 
maximum amount of fat that patients would be likely to consume 
during acute illness and treatment. 

A 26% higher mean absorption rate of piperaquine (95% CI, 
5.1 to 42%) was seen in patients who received study treatment 
with concomitant food as compared to fasting patients when us- 
ing a full covariate approach. This is most likely a result of the 
rapid gastric emptying associated with consumption of liquids 



TABLE 2 Population estimates of the final model describing piperaquine population pharmacokinetics in fasting {n ■ 
with uncomplicated P. falciparum malaria in Thailand 

Value 



15) and fed [n = 15) patients 



Parameter 



Population estimate" 
(% RSE'') 



95% CI for population 
estimate'' 



IIV (% CV)" 
(% RSE'') 



95% CI for IIV*" 



Typical parameters 
CLIP (liters/h) 
VJF (liters) 
Qi/P (liters/h) 
VpJF (liters) 
QJF (liters/h) 
Vp^lF (liters) 
MTT (h) 

No. of transit comp. 

F (%) 
a (% CV) 

Covariate effects 
Dose effect on f(%) 
Age effect on Vpj (%) 



67.6 (11.6) 
3,030 (16.4) 
408 (15.0) 
6,240 (14.6) 
109 (13.6) 
24,400 (10.1) 
2.04 (7.50) 

3 (fixed) 
100 (fixed) 

30.7 (4.42) 



25.3 (34.4) 
4.10 (18.1) 



54.0-85.5 

2,160-4,180 

309-557 

4,890-8,530 

83.3-143 

20,000-29,500 

1.80-2.41 



27.8-33.5 



24.4 (26.0) 
51.6 (32.3) 

45.6 (48.8) 
25.8 (48.2) 

24.1 (52.7)/39.4 (22.7)' 
48.8 (16.6)' 



9.82-53.2 
2.38-5.32 



17.4-29.6 
31.2-68.1 



18.8-68.4 
6.67-37.9 



8.77-35.6/29.2-48.0 
38.3-56.0 



Computed population mean values from NONMEM. Interindividual variability (IIV), between-occasion variability, and random residual variability are calculated as 100 X 
-^£XP(o-^)-l. 

'' Assessed by the nonparametric bootstrap method (n = 1,000 iterations) for the final pharmacokinetic model. Relative standard errors (RSE) are calculated as 100 X (standard 
error/mean). Ninety-five-percent confidence intervals are displayed as the 2.5 to 97.5 percentiles of bootstrap estimates. 
' Between-occasion variability. 

CL, elimination clearance; V^, central volume of distribution; Q, intercompartment clearance; Vp, peripheral volume of distribution; MTT, mean absorption transit time; No. of 
transit comp., number of transit compartments; F, oral bioavailability; ct, additive residual error; CV, coefficient of variation. 



2056 aac.dsm.org 



Antimicrobial Agents and Chemotherapy 



Food Effects on Piperaquine Pharmacokinetics 



TABLE 3 Post hoc estimates of the final model describing piperaquine population pharmacokinetics in fasting (n = 15) and fed (n = 15) patients 
with uncomplicated P. falciparum malaria in Thailand 



Median value for group (IQR)° 



Parameter for post hoc estimate*" 


Total 


Fasting patients 


Fed patients 


P value 


CL/P (liters/h/kg) 


1.02 (0.802-1.21) 


1.09 (0.745-1.31) 


0.988 (0.867-1.10) 


0.548 


V/F (liters/kg) 


516 (378-613) 


488 (375-607) 


517(411-594) 


0.917 


ti,2 (days) 


19.5 (16.8-21.8) 


19.8 (16.2-21.3) 


19.1 (18.1-21.8) 


0.520 


Cmax (ng/ml) 


237 (169-347) 


242 (165-319) 


231 (177-395) 


0.663 


r„ax (h) 


3.57 (3.08-4.27) 


3.36 (2.80-4.20) 


3.64 (3.25-4.16) 


0.520 


AUCi„f(hX |xg/ml) 


26.8 (22.7-34.2) 


23.9 (21.6-36.7) 


27.5 (25.2-32.7) 


0.419 


Day 7 concn (ng/ml) 


29.8 (26.1-35.4) 


26.9 (22.1-42.7) 


30.9 (28.3-34.3) 


0.395 


" Post hoc estimates were calculated as 


median values with interquartile rang 


es (IQR) fi-om empirical Bayes estimates, and statistical differences were estimated with 


I nonparametric 



Mann-Whitney test. 

^ CL, elimination clearance; F, oral bioavailability; V, apparent total volume of distribution { + Vpj + Vp2); fi/2, terminal elimination half-life; AUC^f, area under the 
concentration-time curve from time point 0 to day 138; C^^, predicted peak concentration; T^^, predicted time to peak concentration; Day 7 concn, predicted concentration at 
day 7. 



(28) and/or the increased dissolution of piperaquine when coad- 
ministered with fatty liquids. However, dihydroartemisinin is pri- 
marily responsible for the initial parasite-killing effect, and this 
trend of an increased absorption rate of piperaquine is not likely to 
have any substantial clinical impact with respect to treatment out- 
come. Indeed, the prophylactic protective effect of dihydroarte- 
misinin-piperaquine in volunteers (n = 800) at risk of malaria 
infections in Thailand was not significantly different when admin- 
istered with or without food (29). Furthermore, the trend of an 
increased absorption rate associated with concomitant food in- 
take did not result in higher predicted peak levels of piperaquine 
in the final model, so it is not likely to have any clinical impact. 
This suggests that the manufacturer's recommendation of taking 
the treatment without food (up to 3 h after the dose) could be 
reviewed. 

Evaluation of population versus noncompartmental ap- 
proaches. The modified Monte Carlo mapped power approach 
was implemented successfully and enabled a numerical approach 
to evaluate the power to detect various degrees of food effects with 
the given study design. It demonstrated a marked increase in sta- 
tistical power when using a population approach compared to a 
noncompartmental analysis with a traditional groupwise statisti- 
cal comparison (35% versus 81% difference needed for statistical 
significance with the competing approaches). This evaluation also 
lends further support to the results, since it demonstrated that the 
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FIG 5 Modified Monte Carlo mapped power approach, showing the statistical 
power to detect various degrees of covariate effects with the final pharmaco- 
kinetic model and the given study design (15 fed and 15 fasting patients). 



study was adequately powered to detect clinically relevant phar- 
macokinetic differences. A population approach also offers other 
advantages, such as a mechanistic understanding of the pharma- 
cokinetics and influential covariates as well as a quantitative mea- 
sure of the variability between and within patients. However, it is 
of course more time-consuming than a traditional noncompart- 
mental approach. Small pharmacokinetic studies are often a real- 
ity in antimalarial drug research, and the presented novel power 
approach could be a useful tool to evaluate the power of detecting 
true covariate effects in mixed-effects modeling of a given study 
design. 

In conclusion, piperaquine population pharmacokinetics were 
well described in fasting and fed patients with uncomplicated P. 
falciparum malaria in Thailand by a three-compartment distribu- 
tion model with a flexible absorption model. The final model 
showed a marked increase in bioavailability during recovery from 
malaria but demonstrated no clinically significant impact of con- 
comitant intake of a low-fat meal. The novel methodology for 
assessing the power of detecting a true covariate relationship pre- 
sented may be used in mixed-effects modeling to provide a nu- 
merical approach to conclude a lack of covariate effects. 
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